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SUMMARY 


A compressible linear stability theory is presented for nonparallel 
three-dimensional boundary-layer flows, taking into account the normal 
velocity component as well as the streamwise and spanwise variations of 
the basic flow. The method of multiple scales is used to account for the 
nonparallelism of the basic flow, and equations are derived for the spatial 
evolution of the disturbance amplitude and wavenumber. The numerical proce- 
dure for obtaining the solution of the nonparallel problem is outlined. 


I . INTRODUCTION 


For laminar-flow vehicles, a high performance is achieved with 
respect to range and economy of vehcile operation by reducing the friction 
drag. The design of such vehicles is strongly influenced by the stability 
considerations of the boundary layer. To maintain laminar flow with 
minimum external power, an optimum amount of suction, or cooling (in air), 
or heating (in water) is required. For the design of swept LFC wings of 
transonic aircrafts, this optimization process needs accurate computations 
of the stability characteristics of three-dimensional, compressible 
boundary-layer flows. 

An extensive treatment of the stability theory for compressible 
flows is given by Mack (1969) for two-dimensional mean flows, where the 
disturbance can be two or three-dimensional. These stability theories 
treat the mean flows as quasi-parallel. Some incomplete attempts to 
account for the nonparallel flow effects by including either the normal or 
some of the streamwise derivatives of the mean flow were given by Brown 
(1967), Gunness (1968), and Boehman (1971). 

For two-dimensional heated boundary-layer flows, El-Hady (1978) 
and El-Hady and Nayfeh (1978) introduced a complete nonparallel stability 
theory to account for the rate of heat transfer between the fluid and the 
wall. The nonparallel stability results are in better agreement with the 
heated water experimental data of Strazisar et al (1977) and Strazisar and 
Reshotko (1978), than the parallel results of Lowell (1974). 
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Recently, El-Hady and Nayfeh (1979) analyzed the effect of the non- 
parallelism of the mean flow on the stability characteristics for two- 
dimensional subsonic and supersonic flows. Results calculated by the non- 
parallel stability theory are in better agreement with the supersonic 
experimental data of Laufer and Vrebalovich (1960) and Kendall (1975) than 
the results calculated by the parallel theory of Mack (1969). 

The propagation of three-dimensional disturbances in three-dimensional 
compressible boundary layers was numerically investigated by Mack (1979) and 
Lekoudis (1979). Their analysis was for parallel flows over an infinite 
sweptback wing. Their results show that the effects of compressibility are 
negligible near the leading and trailing edges (regions of cross-flow type 
instability). However, away from the leading and trailing edges (regions 
of T-S type instability), the maximum amplification rate is reduced and 
the most unstable-wave orientation is considerably changed due to compres- 
sibility effects. 

In this article, a compressible linear stability theory is presented 
for three-dimensional disturbance in a nonparallel three-dimensional 
boundary-layer flows. Section II contains the formulation of the problem. 
Section III contains the method of solution for the zeroth and first-order 
problems. The computational procedures are outlined in Section IV. 
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II. FORMULATION OF THE STABILITY THEORY 


We consider the spatial, three-dimensional stability of laminar 
compressible three-dimensional steady viscous flows to small-amplitude 
disturbances . 

The flow field is described by the Navier-Stokes , energy, and state 
equations. Lengths, velocities, and time are made dimensionless using a 
suitable reference length L*, the freestream velocity U* and L*/U*, 

00 CO 

respectively. The pressure is made dimensionless using p*U*^ . The 
temperature, density, specific heats, viscosity, and thermal conductiv- 
ity are made dimensionless using their corresponding freestream values. 

2 . 1 Disturbance Equations 


To study the linear stability of a steady three-dimensional, bound- 
ary-layer flow (basic flow), we superpose a small time dependent disturb- 
ance on each mean-flow, thermodynamic, and transport quantity. Thus, we 
let 


q(x,y,z,t) = Q (x,y,z) + q(x,y,z,t) 


( 1 ) 


where Q^(x,y,z) is a three-dimensional basic-state quantity and q(x,y,z,t) 
is a three-dimensional unsteady disturbance quantity. Here, q stands for 
the velocity components (u, v, and w) , temperature T, pressure p, density 
p, and viscosity y . Substituting (1) into equations governing the flow 
field, subtracting the basic-state quantities, and linearizing the result- 
ing equations in the q's, we obtain the following disturbance equations: 


(p u + pU ) + |—(pv + pV ) + |— (p w + pW ) = 0 
3t 3x s s Sy s s 3z s s 


( 2 ) 
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(6) 


and the state equation. Here R = ^e^e^^*^e Reynolds number, 

Pr = C V /k is the Prandtl number, Ec = U^/C„ T is the Eckert number. 
Pee e P e 

e e 

and (j) , the perturbation dissipation function, is defined as 


3U - 3V 

2r(^ |ii + 


c „ 3U . „ 

S dV ^ S 3 Wn , „ r S fOV dWv 

+ + 2m[^ + _) 


3V 
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(7) 


Moreover, r and m are given by 


r = I (£ + 2 ), m = I (£ - 1) (8) 

where £ is the ratio of the second to the first viscosity coefficients. 

Equations (2)-(6) represent the stability equations for a three- 
dimensional disturbance in a three-dimensional basic flow. 
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2 . 2 Boundary Conditions 


The analysis presented here is applicable to cases with suction or 
blowing as well as cases with uniform or non-uniform wall heating or 
cooling. The stability problem is completed by the specification of the 
appropriate initial and boundary conditions. We consider next the bound- 
ary conditions. 

At the wall, we require the vanishing of the component of the relative 
velocity of the fluid/solid surface that is parallel to the interface, even 
in the region of perforations. This is a reasonable assumption provided 
that the percentage of the permeable area is small and most of the flow 
there is directed normal to the wall. The normal velocity and thermal 
boundary conditions at the wall need careful consideration. 

For an impermeable wall, both the mean and disturbance velocities 
normal to the surface must vanish, whereas for a permeable wall this is 
not the case. The mean normal velocity component is not zero. For dif- 
ferent disturbance normal velocities at the wall, Gaponov (1971, 1975) 
showed a destabilizing effect of a nonzero value for the normal component 
of the disturbance velocity at the wall. His results are based on the 
calculation of the neutral stability curves. Such curves are not directly 
related to transition. Moreover, the boundary condition used by Gaponov 
is not of a practical application. Of more interest are the results of 
Lekoudis (1978), who examined the effect of the normal component of the 
disturbance velocity on the growth rates of disturbances. This boundary 
condition is given in the form of an admittance calculated for different 
configurations. He concluded that the condition of zero normal velocity 
at the wall is a reasonable approximation when the surface permeability 
is very small. 
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The thermal boundary condition for the disturbances needs an analysis 
of the heat conduction problem in the region very close to the wall (e.g., 
Dunn and Lin, 1955). This analysis results in a thermal boundary condition 
for the disturbance in the form 

(0) + g T(0) = 0 

3y 

where g is a constant that depends on the disturbance frequency, and the 
physical properties of the liquid and the adjacent solid wall. This con- 
dition holds for very low frequencies because the thermal fluctuations can 
penetrate large distances into the solid wall. On the other hand, for 
very high frequencies, the thermal inertia of the solid makes the thermal 
fluctuations die out in the solid very close to the surface, and the wall 
remains at the temperature of the basic flow. In this case, the condition 
T(0) = 0 is a very accurate approximation. 

In the freestream, it is assumed that all disturbance quantities die 
out for subsonic disturbances and satisfy a radiation condition for super- 
sonic disturbances. For boundary layers in subsonic flows, disturbances 
have amplitudes that decay exponentially in the freestream. For boundary 
layers in supersonic flows, we restrict our analysis to subsonic disturb- 
ances, that is to disturbances that move subsonically with respect to the 
freestream. The amplitudes of these disturbances decay also exponentially 
in the freestream. 

In this analysis, we consider walls of small permeability and subsonic 
disturbances of sufficiently high frequencies. With these assumptions, the 
disturbance boundary conditions become 

u=v=w=T=0aty=0 (9) 

u,v,w,T -> 0 as y (10) 
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2 . 3 Nonparallel-Flow Considerations 


We consider weakly nonparallel flows. Thus, to account for the non- 
parallelism of the mean flow, we require that all mean-flow variables be 
weak functions of the streamwise and spanwise positions. Moreover, we 
require that the normal velocity component be small compared with the other 
velocity components. These assumptions are expressed mathematically by 
writing the mean-flow variables in the form 

Ug = Ug(xj,y,Zj), Vg = eVg(xj ,y,Zj) , Wg = W^(xj,y,Zj), 

Pg " " T^(xj,y,z,) (11) 

where 

Xj=ex,Zj=Ez (12) 

and e is a small dimensionless parameter characterizing the nonparallelism 
of the mean flow. In the parallel-flow approximation, Vg = 0 and all 
variables in Eqs . (11) are independent of Xj and Zj. The main idea behind 
the nonparallel-flow analysis is to make a perturbation about the parallel- 
flow solution (Bouthier, 1973; Caster, 1974; Nayfeh et al, 1974). The method of 
multiple scales is used to effect this perturbation. Accordingly, 
different streamwise scales x, x^, z, and Zj are introduced. The fast 
scales X and z are used to describe the relatively rapid streamwise and 
spanwise variations of the traveling wave disturbance. The slow scales 
Xj and Zj are used to describe the relatively slow variations of the mean- 
flow quantities', the disturbance wavenumber, the growth rate, and the 
amplitude . 
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III. METHOD OF SOLUTION 


To determine an approximate solution to Eqs. (2)-(ll), we use the 
method of multiple scales and seek a first-order expansion for the disturb- 
ance variables u, v, w, p, and T in the form of a traveling harmonic wave, 
that is 

•1 ( ^ 1 j y j ^ I ) t ; E ) — ^ 1 j y > 2 j ) + . . . jexp(i9) (13) 

where the phase function 0 is defined by 

Ve = K() = [a(,(xi,zi), 3o(xi,zi)J (14a) 

30 

3^ = - (14b) 

Here and Bj, are the quasi-parallel components of the wavenumber vector 
Kq in the x and z directions and 0) is the circular frequency. 

For the case of quasi-parallel spatial stability 0 ) is real and aQ 

and Bo are complex, while for the case of quasi-parallel temporal stability, 
and Bp are real and oj is complex. For the general case, a^, Bj, and co 
are complex. The present study is limited to spatial stability so that W 
is a known real quantity. 

The viscosity disturbance is related to the temperature disturbance by 

,y,Zj ,t;e) = -^ T(xj ,y,zi ,t;e) (15) 

s 

while the density disturbance is related to the temperature disturbance 
by the state equation. 
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Substituting Eqs. (13)-(15) into Eqs . (2)-(6), transforming the time and 


spatial derivatives from t, x, and z to 6 , Xj , Z| , and equating the 
coefficients of e° and e on both sides, we obtain ' problems describing the 
q^ and q^ disturbance quantities. These problems are referred to as the 
zeroth and first-order problems. They are given next for the physical 
problem of air boundary layer treated as perfect gas. 


3 . 1 The Zeroth-Order Problem 

Substituting (13)-(15) into (2)-(ll) and equating the coefficients 
of on both sides, we obtain 

h j ( u u , V Q , w Q , p Q , Tq ) — ip^CcHgUij + B|jWg) + i(ct(jU^ + Po^g ~ 

+|—(pvg)=0 (16) 

dy s “ 
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1 


3Ug 

i f 

3Vg 

1 


dv 

s 
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1 

3^Ug 

r 

9y 

3y 

R 

3y 

r 

dT 

s 

9y 

3y 

~ R '^s 

3y^ 
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( 21 ) 


^ 0 as y 


( 22 ) 



The operators correspond to the continuity, x-momentum, 

y-momentum, z-momentum, and energy equations, respectively. 

Equations (16)-(22) constitute an eigenvalue problem that is 
solved numerically. The numerical solution is discussed in the next 
section, and the solution of Eqs. (16)-(22) can be expressed in the form 


Uo 

= A(xi , Zj ) 

Cl 

(xj ,y,zi ) 

(23a) 

Vq 

= A(xi ,Zi ) 

C 3 

(xi ,y ,zi ) 

(23b) 

Po 

= A(xj ,Zj ) 


(xj ,y,zj ) 

(23c) 

To 

= A(xj ,Zj ) 

^5 

(Xj ,y,Zj ) 

(23d) 

^0 

= A(xj , Zj ) 


(xj ,y ,zj ) 

(23e) 


The amplitude function A(x.j ,z^ ) is determined by imposing the solvability 
condition at the next level of approximation. 


3.2 The First-Order Problem 


Substituting (13)-(15) into Eqs. (2)-(ll) and equating the coefficients 
of £ on both sides, we obtain 


L.(u,,Vj 

,Wj ,Pj ,Tj ) 

= I. 
1 

for i = 1 , 2 , . . , 5 

(24) 

U^ = Vj 

= -1 = T, 

= 0 

at y = 0 

(25) 

Uj ,Vj ,Wi 



as y -*■ 00 

(26) 


where the operators L 1 -L 5 are defined by Eqs. (16)-(20) and Ii-Is are the 
inhomogeneous terms in the continuity, x-momentum, z-momentum, and energy 
equations, respectively. These inhomogeneous terms reflect the effects of 
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the streamwise and spanwise variations of the disturbance amplitude, the 
normal basic-velocity component and the streamwise and spanwise variations 
of the wavenumber. Using Eq. (23), these inhomogeneous terms are defined 
by 
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f edx.^ R Os dx^ dy dz^ R Osdz^ dx^ 

1 3U 3C 

} C, + f y (Y - l)M^ f-r T — + 

1 R 55 p dv 3x 


T 3x, 
s 1 


R s 

y2 


3T 


3W 3? YM 3T 

+ m T— ] - ~ (U „ ... 

3z, 3y T s 3x s 3y 
1 s 1 

9^„ 


^ ■ V W 


5”, 

3U 

(m -r— ^ 

3x, 

9V 

3y 3Zj 


9T 

3z > ^ 

+ (Y - 

1)M^(U 

e 


3C 


3^ 3T 3T 3T . 3a 

+ V -r— + W ^ ) + [;t 2 - ■ (U -r—- + V ^ + W '^— ) + (p ^ 

s 3y s 3z T s 3x, s 3y s 3z, RPr s 3x, 

•' 1 s 1 1 1 


+ 2a 


0 3x 


~ (u 


33. 


3y, 


RPr s 3z 


+ 20 


)] c, + ( 


2 i 


0 3z " ”5 
1 




S ) 


RPr o^s Ts 3x 


W 9? V 3i;^ 
, ^ 2l S. 5 s 5 

^RPr T ^ 9z ' T 3y 

s 1 s 


9 P 

{<v - ‘KCar 


1 

9U 


a u ( 


o^s '9z 


3W 3U 3V 3W 

Ss 2 i „ . s s , s 

+ T- — ) + -— 3 „y (m -5 + m -5 + r — 

dx^ R Os dx^ dy dz 


1 s 1 ; 


(31) 


The inhomogeneous Eqs . (24)-(26) have a solution if, and only if, a 
solvability condition is satisfied. This condition demands the inhomo- 
geneous terms be orthogonal to every solution of the adjoint homogeneous 
problem. In the next section, we obtain equations for the modulation of 
the wave amplitude and the wavenumber with position using the solvability 
condition . 
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IV . COMPUTATIONAL PROCEDURE 


The zeroth-order stability problem defined by Eqs. (16)-(22) is 
reduced to the solution of the following set of eight first-order differen- 
tial equations with eight homogeneous boundary conditions 


• 8 


where 


Ll 

8y 

- Z 
j = 

a . . < 

1 


for 

1— 1 
It 

2, .. 

8 


(32) 

Zol = 

= Zo 3 

- Zq I 

5 = Zo7 = 

= 0 

at y = 

0 



(33) 

Zo 1 ) 

Zo 3 > 

Zq 5 > 

Zo 7 ^ 0 

as 

y CO 




(34) 

Zo 1 • 

= Oo , 

Zo 2 

9uo 
3y ’ 

Zo 3 

= Vo , 

Zo 4 = 

Po , 



Zo 5 = 

= To 

Zo s 

3To 

- 8y * 

Zo 7 

= Wo , 

Zo 8 = 

9wq 

9y 


(35) 

le a . 

. are 

the ( 

Blements 

of 

8x8 var 

iable- 

coefficient matrix. 

The 



ij 


nonzero elements of this matrix are 


ai 2 = 1 





321 _ iR 

(aoU 

s 

+ 3oW 

s 

- cj) 

+ ao + 

T y 
s s 

1 

^2 2 ” 

9y 




R 

9U 

s 

ia.t-J- 


. f 

32 3 .j, 

s s 

9y 

9y 

T 9y 

S 


aai* - fyM^aoCaoU + 3oW - co) 

u ess 

s 


2 

0 
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fC6o 


az 5 ” 


(ctoU + 3oW “ <Ai) 
s s 


1 9 


du 9U 

(—1 _±) 

Ug 9y dTg 9y 


32 6 


V, 3y 


33 1 ~ 


33 3 “ 


1 

1 s 

T 9y 
s 


334 = " iyM^(aoU + 3 oW - co) 

e s s 


335 = ^ (aoU^ + 3oW^ 

s 


&i) 


337 

34 1 


iBo 


9 y 


9T 


/ ^ s r s X 
IXO.0 9y ^ 


34 2 

= - iycto 












9y 

9T 

9^T 

iR 

to)] 

34 3 

= X [-cto - 

g? + 

r 

T y 

S £ 

s 

, 3y 
} 



9y 

r s 

^ T ^ 
s 

--T ta.u, -BoW, - 

S S 




By 

T 

9T 


9U 

9W 

344 

= - ixryM^ 

[(— 

s 

9y 

1 

+ — *• 

T 

s 

1) 

9y 

(aoU + 
s 

3 qW^ - (ii) + Qo Qy ■*■ Bo 




dy 



9U 

9W 


- to)] 

34 5 

s 

s 

dT 

s 

s 

(ao 


^0 

+ ^ (ctoU + BoW 

T y B y s s 

s s 


346 - (CtoU^ + 3oW^ - 0)) 

s 

9y 9T 

a47 = - ix3o(— ^ + f- 

S s 
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I 


34 8 
356 

36 2 
36 3 

3e 4 

36 5 

30 6 

36 8 
38 3 

3s 4 

38 5 

38 6 

33 7 


= - ixBo 


= 1 


3U 

- 2Pr(Y - DM^ ^ 


RPr 

Vs 'y 


au 


aw 


2iPr(y - DM^Oot;:— + 3o 


e oy 


3y 


(Y - DM^ (aoU + 60 W - (o) 

y ess 


iRPr 

T y 

s s 


dy au 

(aoU^ + 3oWg - (o) + ao^ + Bo" - (y-Dm^ ~ ~ 

^ s s ^ 


aw , a^y 

_ i i 

^3y y ay2 

s 


2 




aw 

- 2Pr(Y - DM^ -V 
e 3y 


.Q /1 s r s., R s 

“ V ^ ^ ^ V~T~ 

s s s s 


. iRBo 


fBoYM^ (aoU + BoW - co) 

6 S S 


efl 1 c. 

" ®«“s - “> - 3^ <sr af > 

S S S 


, dy aw 

1 s s 


y dX 8y 
s s 


(otoU + BoW - Cj) + ao + Bo 
i y s s 
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1 

-BB - - y, ay 

where 

X = l/[^ + iryM^CaoU^ + BqW^ - co)] 
s 

Equations (32)-(34) constitute an eigenvalue problem and it has nontrivial 
solutions only for certain combinations of the parameters , 3o » ^ and R. 


4 . 1 Eigenvalues and Vectors 


Outside the boundary layer (at y = y^ 
layer thickness), the mean-flow quantities 
nonzero elements of the coefficient matrix 
constants. They become 


> 6 , where 6 is the boundary- 
are independent of y and the 
^a^ j J given by Eq. (32) are 


ai 2 “ 1 ) 32 1 “ iR(cto + CBo “ co) + Cto + Bo > 

32 4 = iROo - fYMgOto (Olo + CBo “ . 

32 5 = fcto (ao + CBo - cu), 33 1 = - iao , 334 = - iyM^(cxo + CBo “ cd) , 

S3 5 = i(ao + CBo - co), S3 7 = - iBo , S42 = - iX“o , 

S43 = - X [iR(cto + CBo - co) + ao + Bo ] , aoo = ixr(ao + CBo “ co), 
aoe = - ixBo, ^ss = 1 , as 4 = " iRPr(y - l)M^(ao + CBo - co), 

Ses = iRPr (ao + CBo ~ co) + ao + Bo j as 4 = iRBo ~ fyM^Bo (c^o + CBo ~ co). 


as5 = fBo (oto + CBo "■ co) , Ss? - a 2 i 


20 



where 


X = 1/[r + iryM^Cao + CBp ~ <^)3 » and C = at y = 

Equations (32) with a constant coefficient matrix [a^^J permit a solution 
that can be expressed in the general form 


8 


Z A . . c . exp(X ^y ) for i = 1,2,..., 8 


Zo . 

" j=l J 


(36) 


where the X. are eigenvalues of the matrix fa..!, the A .. are the 
J ij ij 

elements of the corresponding eight eigenvectors, and the are arbi- 
trary constants. 

The values A. and A. . can be derived analytically by rewriting the 
J ij 

eight first-order equations (32) with constant coefficients as four second- 
order equations in the form 


d^J 

-■■■■ 7 ' - Z b..J. = 0 for i = 1,2, ..,4 

dy j=i j 


(37) 


where 


Jl - Zo 1 , Jz - Zo 4 , J 3 ~ Zo 5 > J4 - Z 


0 7 


and the coefficients b. . are given by 


b 

1 1 

= 321, b 

1 2 

~ 32 4 > bi 3 - 

A /S A /N 

32 5 , 

A A 

ba 2 

- a24342 

/s 

+ 

334343 + a463G4 

A A A A 

+ 

34 8 3a 4 

A A 

ba 3 

= 325342 

+ 

335343 + a46ag5 

+ 

348 3a 5 

bs 2 

=^64. b 

3 3 

= 365, bo 2 = 

38 4 J b4 3 


(38) 
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The solution of Eqs. (37) has the form 


J. = E B . .d ,exp(X .y) for i = 1,2,..., 4 
1 ij J J 

where the are the same as the eigenvalues of Eqs. (32) with constant 
coefficients, the are the elements of the corresponding eight eigen- 

vectors, and the d^ are arbitrary constants. From the characteristic 
determinant, it follows that the eigenvalues are 

Xi,5 = + (bii) ^ ^ 


X 2 ) 6 - + 


I'J (b22 + bss) -t- (b22 “ t>33)^ + b23b32] | ^ ‘ 

,7 = + I'J (b22 + bss) - (b22 “ t>33)^ -t- b23b32] | 


V 2 \ V 2 


Xif J 8 - + (bi 1 ) 


V 2 


The B.. can be obtained from the solution of the characteristic 
ij 

equation. They are given by 


B.. “* 1, B 2 . “ O 3 ^> 4 . — 0 


for j = 1,5 


a 2 ^ \ /N /\ /N 

. “ 3651324 + 3253GJ 




(321 - X^) 
J 


, B 2 . - 36 5 “ X j , 


B 3 . 
J 


36 4 3 B4 j 


364385 + (X. - 355)384 

(321 - X^) 


for j = 2,3,6, 7 


( 39 ) 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 


22 



and 


Bij = 0, B2j = 0, Ba^ = 0, B4^ = 1 


for j =4,8 


( 46 ) 


The A. . are related to the B. . by Eqs . (36)-(39): they are 
ij ij 


Aj. - 1 , Aa - - A., A3. - (aaiBi. + aa^Ba. + S35B3. + 33764. )/A.Bi., 

J JJJ J J J JJJ 


A 4 . - Ba./Bi., As. - B 3 ./B 1 ., As. - X.B 3 ./B 1 ., Ay. - B 4 ./B 1 ., 
J JJ1 JJJJJJJ JJ 

As . = (33462 . + assBa . + 39764 .)/X.Bi . 

J J J J J J 


(47) 


These eigenvectors are normalized such that 


Zq 1 = Y. c.exp(X.y) at y = y 
j = l J J 


4 . 2 Boundary Conditions 

The boundary conditions at infinity (34) demand the constants 

C 5 , cs , cy , and ce be zero. To set up these boundary conditions for numerical 

solution, we first solve Eqs. (36) for the c.exp(A.y) and obtain 

J J 


c.exp(A.y) = E f..Zo- for j = 1,2 , ...8 

J J i=l ^ 


(48) 


where the matrix [f ] is the inverse of Setting C 5 =ce=C 7 =C 8=0 in 


Eq. (48) leads to 


Z e. .Zn . = 0 for i = 5, 6 , 7, and 8 at y = y 
i=l 1 "'e 


(49) 


Equations (49) replace the boundary conditions Eqs. (34). 


23 



The boundary conditions at the wall (33) can be set up for numerical 


solution by writing them in the form 


Z e. ,Zi . = 0 for i = 1,2, . . .8 (50) 

j=l ^ 

where the e. . are the elements of an 8 x 8 matrix with only four nonzero 
ij 

elements . 


4.3 Integration and Orthonormalization 


For the spatial stability problem, we assign values to Ci) and R and 

two relations among CXq^, chq ^ , 6oj.j and 3o^, where ao = cto ^ + icto and 

3o ~ + i3o^- Then, we guess the remaining two relations. We determine 

the f. . in Eqs. (49) and use this boundary condition to construct a linear 
ij 

combination of the general solution given by Eqs. (36). As y the four 

growing solutions in Eqs. (36) are eliminated. A variable step size algorithm 
developed by Scott and Watts (1977), based on the Runge-Kutta-Fehlburg 
fifth-order formulas, is used to integrate Eqs. (32) from y = y^ to the wall. 

A straightforward integration fails to produce four linearly independent 
solutions because of the buildup of parasitic errors among the different 
solutions. To overcome this difficulty, the integrator used is coupled with 
an orthonormalization test that is based on the modified Gram-Schmidt proce- 
dure . 

Since testing for independence after each integration step is 
expensive, we use a modified algorithm (Darlow et al , 1977) and choose a 
preselected set of points where orthonormalization is performed. These 
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points are assigned a priori by using information about the points where 
orthonormalization is needed. 

At the wall, the values of the linearly independent solution vectors 
are linearly combined to satisfy all but one of the wall boundary conditions. 
The last wall boundary condition can only be satisfied by this combined solu- 
tion when the exact remaining relations among Cto ^ , 3oj.» and 3o^ have 

been found. A Newton-Raphson procedure is used to determine these relations. 
With the eigenvalue determined to within the desired accuracy, the eigen- 
functions can be recovered using the stored solution vectors. This solution 
can be expressed in the form 

Zo ^ = A(xi,zi) (xi , y,zi) for i = 1,2, ...,8 (51) 


4.4 Solvabil ity Condition and the Adjoint 

With the solution of the zeroth-order problem given by (51), the 
first-order problem becomes 


3Zj 

i s 

— - y 

a Zi - r + F 

3A 

3y 

j=i 

ij j 1 3xi 1 

3zi 

Zi 1 

~ ^13 

= Zi 5 = Zi 7 = 0 at y 

= 0 

Zi 1 

9 3 J 

Zi 5 . Zi 7 -> 0 as y ^ 

oo 


(52) 

(53) 

(54) 


where the , E^, and are known functions of the Uq , 3o » and mean- 

flow quantities; they are defined by 


^ 3A ^ 3A „ 


(55a) 
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(55b) 


Gz 


8A ^ 3A . . ^ R 

■ 5 — + E, - 5 — + DzA = - ifaoT I, - — I, 

9x1 2 3zj ” s 1 ys 


9A „ 8A „ . _ 

_ + E3 ^ + D3A = T^Ii 


(55c) 


^ 9A ^ 3A ^ . , 


'"^.2^)1, H-rxT .Mx3(53d) 

dy 1 <?3vl 11 ^ 


s 3y 


Gs + Es + DsA = 0 

dxi dzi 


(55e) 


Ge 


9A 


9xi 9zi 


Eg + DgA — 


RPr 


Is 


(55f) 


„ 9A 9A 

G 7 ~R — + E? - 5 — + D 7 A = 0 

dxi dzi 


9a 


3A 


3X1 - Ee ^ + DeA = - ifBoT^Ii 




(55g) 

(55h) 


where Ii - I 5 are defined by Eqs . (27)-(30). 

Since the homogeneous parts of (52)-(54) are the same as (32)-(34) and 
since the latter have a nontrivial solution, the inhomogeneous Eqs. (52)-(54) 
have a solution if, and only if, a solvability condition is satisfied. In 
this case, the solvability condition is 


00 8 

/ Z 

0 i=l 

where the W. (x 
1 

corresponding 



9A 

3xi 


+ E 


I ,y,Zj ) are 
to the same 


3A T " 

. T + D.A I W.dy = 0 

1 3 zi 1-1 1 ■' 

solutions of the adjoint homogeneous problem 
eigenvalue. Thus, they are solutions of 


(56) 


3W. 
1 

9y 


+ 


8 

E 

j=l 


a. .W. 
Ji J 


0 for i = 1 , 2 , ... ,8 


(57) 
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(58) 


>N /S 


W 2 = Wi* = We = Wa = 0 at y = 0 
^ 

W 2 j W 4 , Wa , We 0 as y 


(59) 


We solve the adjoint problem, Eqs. (57)-(59), following the same 
numerical procedure used to solve the zeroth-order problem. Outside the 
boundary layer (at y = i Eqs. (57), has constant coefficients and 
its solution can be written in the form 


W. = E A* . c*exp(X .y) for i = 1,2,. ..,8 
1 J J 


(60) 


where the A. are the same as those for the zeroth-order homogeneous problem, 
J 

Eqs. (32)-(34), but the A*, are different from the A... The A*. can be 

ij ij 

obtained analvticallv in the same way we obtained the A. .. The A*. 

3 J ij ij 

components are given by 


At, 


A. . 


where 


1 , Aaj = - (^ij 34 2 B 2 j ) /X jBi . , Asj - - 34 3 B 2 j /X jBi j , 


B2./B1., As . = B3./B1., As- = - (a 4 sB 2 . + B3 )/X.Bi., 
Ill 111 .. jJJ 


B4./B1., As. - ~ (a4sB2 . + B4.)/X.Bi. 
Ill J J J J 


Bi . - 1 , B 2 . 


(X^-aa 5 ) 32 4 +32 5 36 4 


J (a65-X^.)(b22-X^)-aG4b 


2 3 


(61) 




a24b23“(b22“X^)a25 


B 4 . = 0 


(a65“Xj)(b22“Xj)-as4b23 


for j = 1,5 


(62) 
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Bij = 0, B 2 j = -364, BSj - b22“^j, 84 ^ - 0 

for j = 2, 3, 6, 7 (63) 

and 

a 2 ^ \ yN /s /s 

.-36 5 7384 + 385364 

Bi . = 0, B2 . = J— — 

^ ^ ( 36 5 j ) (b2 2~A j )~a6 4b2 3 

as 4b2 3 “(t)2 2~)^^ )aa 5 j. 

B 3 . = ^ , B^. = 1 

^ ( 36 5“^ j ) (b2 2“^ j )~36 4b2 3 

for j = 4,8 (64) 

In solving the adjoint problem, the eigenvalue relations we found 
before are used in one integration to produce the adjoint solution. The 
solution of the adjoint problem can provide an independent check 

on the eigenvalues obtained earlier. Moreover, solving the adjoint problem 
provides an easier and accurate way of calculating the group velocity 
instead of the approximate and lengthy finite difference techniques. 


4 . 5 Amplitude and Wavenumber Equations 


Substituting for G., E., and D. from Eqs . (55) in the solvability 
'=111 

condition (56), we obtain the following equation for the modulation of the 
wave amplitude A with position 


Qi 



Q2 


9A- 

8zi 


Hi A 


where 

00 

Qi = / 

0 


8 

Z G.w.dy 

i=i " " 


(65) 


( 66 ) 
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( 67 ) 


<vj 8 

Q, = / E E.W.dy 

0 i=i ^ ^ 

oo 8 

Hi = / E D.W.dy 

0 i=i 1 ^ 

Here, Qi and Q 2 are proportional to the components of the group velocity 

, 9oj 3ca .. 

^3oto ’ a3o 

To determine Hi, we need to evaluate 9ao/9xi, 3ao/3zi, 33o/9xi, 33o/9zi, 
3^./3xi, and 3c./3zi. To accomplish this, we replace by in (32)- 
(34), differentiate the result with respect to x, , and obtain 


3 3^i 

— ( ^) 

3y ^3xi 


8 

z 

j=l 


a . . 

ij 


'3lf> 


3ao ^3o 

G. + E. + S . for i 

1 dXi 1 dX] XI 


1,2, .,8 (69) 


|?i = |?A = |i^ = |Si=Oaty=0 

3xi oxi dxi dxi 


(70) 


9 ? 1 34 3 34 s 347 

9xi ’ 9xi ’ 9xi ’ 9xi 


as y ^ 00 


(71) 


Similarly, differentiation of (32)-(34) with respect to zi yields 


9y 


34i 


8 

E 

j = l 


34. 

a . . ( ' ^— ^ ) 
ij 9zi 


3ao 
i 3zi 


+ 


E. 

1 


+ S . 
9zi zi 


for i 


1,2,.. .,8 (72) 


34 1 _ 9^3 _ 9g 5 _ 34? _ 
3zi 3zi 3zi 3zi 


at y = 0 


(73) 


34i 343 34s 34 ? 

3zi ’ 3zi ’ 3zi ’ 3zi 


0 as 


y ^ OO 


(74) 
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Here G. , E^, and are known functions of cto > 3o > and the basic-flow 
quantities; they are given by 


8 9a . . 

G. = E C. 


i ‘’j ao ' 3o , XI 9ao for i = 1,2, ..,8 (75) 


3a. . 

U 


8 3a- . 

E Z — ^ 

j = l 


j 38o OLo ,xi 


9Bo 'cio.zi for i = 1,2 , ..,8 (76) 


3a . . 
n 


8 da . . 

JJ. 


S . = E C- 


3a . . 

±l 


^xi 3x1 'cto,3o, ""zi '’j 3zi ‘an ,3o 


for i = 1,2,.. ,8 (77) 


Again, applying the solvability conditions to (69)-(71) and (72)-(74), 
we obtain equations for the modulation of the wavenumber with position 


3ao 

9x1 


9ao 
3zi 


(78) 


36o „ 33o _ „ 

37T ~ 


3z] 


(79) 


where Qi and Qa are given in ( 66 ) and (67) and Ha and H 3 are given by 

00 8 

Ha = i / E S .W.dy 
0 j=l J 


(80) 


H 3 = i / E S .W.dy 


0 j=l 


zj J 


(81) 


The quantities Hi, Ha, H 3 , Qi and Qa in Eqs . (65), (78), and (79) are 
slowly varying functions of x and z. For a parallel mean flow, the H's 
vanish and the Q's are constant., Nayfeh and Padhye (1979) derived equations 
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similar to (65), (78) and (79) for incompressible nonparallel three- 
dimensional flows. 

In the spatial theory ao and 3o are complex and 03 is real. We 
define a real wavenumber vector of magnitude ko and direction \p according 
to 

ko — (etorj^or)3 4^ — tan ( 3 o r^ ttoi^) (82) 

and a real spatial amplification-rate vector of magnitude Oo and direction 
according to 

Oo = (oioi, 3oi). ^ = tan'^(3o i/cto i) (83) 

The solution of the eigenvalue problem, gives the complex dispersion 
relation 

01 = oj(ko ,ao ,x,z) (84) 

For fixed 03 , x, and z, there are four real parameters, ko , O, and l|3 . 

Two of them can be determined from the eigenvalue calculation. 

In general, the direction of the wave propagation is different 
from the direction of the wave amplification ij;. The propagation angle 3j3 
can be used as an input parameter, while the question of determining the 
direction of the amplification ijj is still open. Mack (1977) chose the 
direction given by the real part of the group velocity angle to be the 
direction of amplification. He showed that for two-dimensional basic flows, 
the direction of the group velocity deviates by a few degrees from the 
streamwise direction. This deviation decreases as the Mach number increases. 
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Lekoudis (1979) and Runyan and George (1979) chose the direction of 
amplification to be the direction of the local potential flow. An 
amplification direction given by the real ratio of the complex group 
velocities was derived by Nayfeh (1979). For a parallel mean flow, the 
amplification direction is given by the real ratio of 9ao/33o, which was 
derived by Cebeci and Stewartson (1979) and Nayfeh (1979). For a mono- 
cromatic wave generated by a source oscillating at frequency CO at x = 0 
to penetrate large values of x and z, the ray equation 


4^ = = a real quantity (85) 

dx Qi 

defines the direction of the wave amplification for the physical problem 
of real x and z. The wave amplitude and wavenumber will vary along the 
ray as 


' A-0 

exp( 

/Hi dC ) 

(86) 

= / 

H2 

dC 

(87) 

= / 

H3 

d? 

(88) 


Equations (86)-(88) are derived from Eqs. (65), (78), and (79) by using 


dxi „ , dzi „ 

^ and — = Q2 


(89) 


Using Eqs. (14), and (86)-(88) in Eq. (51), we obtain 

Zo i = Ao?^(Xj ,y,Zj ) exp |^iy^(ao + Bo ^ dx-icotj + 0 (e) 


(90) 


where Zq i is related to the disturbance variables by (35), and the 
constant Aq is determined from the initial conditions. It is clear from 
(90) that, in addition to the dependence of the eigensolution on Xi and zi , 
the amplification of the disturbance is a function of the normal distance 
from the wall. 
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